Thermoelectric transport properties of an apparent Fermi liquid: Relation to an 
analytic anomaly in the density of states and application to hole-doped delafossites 
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Through the motivation of the recent discovery of dispersionless regions in the band structure of 
the delafossites, a model density of states of free fermions including a discontinuity as analytical 
anomaly is studied. The resulting temperature dependence of the chemical potential is obtained 
both exactly and by different approximation schemes which are then discussed thoroughly. This 
includes the introduction of an approximation of the polylogarithm difference which is capable of 
accessing a parameter range neither covered by Sommcrfcld expansion nor by Boltzmann approxi- 
mation. It is found that the Fermi temperature and several other temperature scales may be very 
low, giving rise to experimentally observable behaviours differing from the one described by Fermi 
liquid theory. In particular, two kinds of apparent Fermi liquid behaviour emerge at intermediate 
temperatures. This behaviour is related to recent transport data reported for CuCri_ x Mg x 02 [A. 
Maignan et. al, Solid State Commun. 149, 962 (2009)] and CuRhi_ x Mg x 2 [A. Maignan et. al, 
Phys. Rev. B 80, 115103 (2009)] by means of the temperature independent correlation functions 
ratio approximation. In this way an effective density of states as well as the effective charge carrier 
density of these materials are determined. Furthermore, conclusions about the specific heat of the 
latter material are drawn which presents particular effects of the analytical anomaly. 
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I. INTRODUCTION 

Thermoelectric materials provide new opportunities for 
energy savings with possible usage in a vast area of ap- 
plications. From the fundamental point of view their 
optimization remains challenging, but several classes of 
promising materials have been identified. They can be 
manganite&i, titanates^, clathrates^ or skutterudites^, 
to quote a few, but also correlated metals, such as 
Na x Co02— or band insulators, such as In 2 _ x Ge x 02— . In 
the former case their large thermopower was attributed 
to narrow bands, while sharp band edges seem to play an 
important role in the latter case, thereby motivating fo- 
cusing on lower dimensional materials. In this context 
the delafossites^, such as Mg-doped CuCr02— >^ or 
Mg-doped CuRh02«^, show great promise. Even though 
they were not yet thoroughly studied, their dimcnsionlcss 
figure of merit ZT exceeds 0.1 above 1000 Ki^. 

At those large temperatures usually a temperature in- 
dependent Heikes behaviour— ~— is expected. Neverthe- 
less, many materials were synthesized showing a T-lincar 
thermopower— ~—. This behaviour can be explained by 
Fermi liquid theory, but requires a large temperature 
scale. However, recent findings in manganite s 24 > 25 and 
delafossite s 11 i 12 i 26 show slight deviations from Fermi liq- 
uid behaviour which challenge this interpretation. These 
modifications include the observation of non-vanishing 
zero temperature values of the linearly (quadratically) 
extrapolated thermopower (conductance) and a change 
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in the slope between the low and high temperature be- 
haviour which for instance can be seen in Fig. [T] A struc- 
tural phase transition might explain such a crossover, but 
has never been evidenced experimentally. In order to 
study this phenomenon in more detail, it will be referred 
to an apparent Fermi liquid (AFL) in the following to 
account for the differences to an ordinary one. 

Although much effort has been devoted to the determi- 
nation of the electronic structure of the materials where 
the AFL behaviour was observed— £2a2&rZL£L2& little is 
known of a microscopic mechanism creating such a be- 
haviour. However, a recent study— revealed a dispersion- 
less region in the vicinity of the Fermi energy along a par- 
ticular direction of the Brillouin zone. Conclusively, this 
leads to a discontinuity in the density of states. While on 
the one hand such features might be difficult to resolve 
numerically, physical properties on the other are known 
to be strongly dependent on such analytical anomalies. 
As found for other anomalies, like a cusp in the density 
of states^ or van-Hove singularitie s 31 ' 32 , this could ren- 
der other phases observable, like marginal— or non- Fermi 
liquid— ones. Therefore in the presence of a discontin- 
uous density of states novel phases should emerge which 
might explain an AFL behaviour. 

With the purpose of investigating the implications of 
such a discontinuous density of states, in particular an 
AFL behaviour, a phenomenological framework is ad- 
dressed. In this framework the chemical potential can 
be evaluated by the exact solution of the grand canoni- 
cal system using polylogarithm functions^ 3 , and applying 
different approximation schemes to determine the tem- 
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Figure 1: Thermopower S and resistivity a~ 1 (inset) 
of CuRho.9Mgo.1O2 as functions of the temperature T 
(circles). The linear, respectively quadratic, regions are 
highlighted by the dotted lines in order to show the tran- 
sition from a behaviour, which can be explained within 
a Fermi liquid picture, towards a behaviour which is re- 
ferred to as the one of an apparent Fermi liquid (AFL) 
and characterised by the additional offsets So, 5q~ . 



perature dependence of the chemical potential. In order 
to access the region between Sommcrfcld expansion and 
Boltzmann approximation this includes the introduction 
of an additional approximation. While this approxima- 
tion starts from the difference of the polylogarithm func- 
tions, it can be interpreted as a way to obtain an av- 
eraged density of states and a renormalised doping for 
a given temperature even for general density of states. 
Based on the chemical potential of the model density of 
states, the mechanism allows to identify the various tem- 
perature ranges associated to the different behaviours. 
This includes not only the Fermi liquid behaviour, but 
also two types of AFL behaviour with different appar- 
ent effective masses. The evaluation of the thermopower 
will be performed using the temperature independent cor- 
relation functions ratio approximation (TICR) 12 i 25 . It 
enables to access this quantity within the phenomcno- 
logical approach in contrast to the ones starting from 
model Hamiltonians, like the high frequency formulation 
of transport coefficient s 34 ' 35 or the dynamical mean field 
theor y 36 ' 3 ' . Additionally, the method, generalizing mod- 
ified Heikes formula s 18 i 38 ~ ^ can account for some dis- 
order and interaction effects, contrary to approaches us- 
ing the constant scattering time approximation within 
Boltzmann transport theor y 12 ' 21 ~ 23 ' 41 . Furthermore, the 
TICR approximation allows to discuss the thermopower 
when governed by its thermodynamic part, relating it to 
the chemical potential. In order to gain further knowl- 
edge of the new phases, the model is then applied to 
CuCn_ x Mg x 2 -ii and CuRhi_ x Mg x 2 -^. Thereby ex- 
tracting not only the parameters of the density of states 
out of the thermopower, but accessing the effective charge 
carrier density as well. The obtained results are then 
compared to the ones given by the calculation using first 



principle studies and by statistical electron diffraction 
spectroscopy (EDS). 

The structure of this investigation concerning the rela- 
tion of an AFL behaviour and analytic anomalies starts 
with the outline of the general framework in Sect, [ill 
Therein the form of the used effective density of states 
will be motivated for the delafossites in more detail and 
effective variables are defined. In the next part, Sect. IIII1 
theoretical background calculations are given in order to 
determine the phase diagram based on the behaviours 
found for the chemical potential. Therefore it contains 
the discussion of different approximation schemes. On 
the one hand Sect. IIII.AI places Sommerfeld expansion 
in the framework of polylogarithm and exemplifies the 
determination of characteristic temperature scales, while 
on the other hand the results of Boltzmann approxima- 
tion are pointed out in Sect. IIII.BI and are prepared for 
the observation of an AFL behaviour within this regime. 
In addition a third approximation examining the region 
between these two by the means of the approximation 
of the polylogarithmic difference (APLD) is introduced 
therein, too (Sect. IIII.C|> . It is then followed by Sect. |IV] 
including a detailed discussion of the TICR framework 
and the implication to the thermopower of these approx- 
imation schemes fSect. lTVA"| as well as the application of 
the approach to the thermopower data of CuCri_ x Mg x 02 
fSect. lTVBl and CuRhi_ x Mg x 2 (Sect. ITVCj) . The spe- 
cific heat following from the obtained effective density of 
states of the latter material is discussed therein, too. The 
results are summarized in Sect. |Vj 

II. GENERAL OVERVIEW 

In order to find a suitable approach for the delafossites 
showing an AFL behaviour their band structure should 
be considered. First principle studie o 11 ' 21 ' 28 showed a 
sharp increase in the density of states at the Fermi en- 
ergy close to the upper band edge. Since these inves- 
tigations did not focus on the analytical behaviour of 
this quantity, possible discontinuities and their implica- 
tions were not addressed. A more recent examination of 
CuRhi_ x Mg x 02 indeed pointed out non dispersive bands 
in particular directions in the vicinity of the Fermi en- 
ergy—. Therein it was proposed that this could be the 
origin of degenerate states causing the AFL behaviour. 
This motivates investigating a density of states which is 
dominated by a discontinuity po at the band edge. Addi- 
tional terms as appearing in a Taylor series^ of a single 
band 

00 

p{e) = Y,^e n -Q(e)®(W-e), (1) 

n=0 

where W denotes the band width, should therefore be 
small corrections. With such a density of state an AFL 

*In fact, the calculation can easily be generalized to arbitrary 
(non-integer) powers in the scries by replacing the factorial by the 
Gamma function. From this procedure complications only arise 
when taking into account a finite band width W or in the low 
temperature expansion. 
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behaviour should accurately be described since at the ob- 
served large temperatures the chemical potential is be- 
lieved to be located in the band gap. Then the product 
of the Fermi function with the density of states changes 
mostly at the band edge. In the finite temperature cal- 
culations following, mainly a sufficiently large band is as- 
sumed, i.e. W — > co, while effects of a finite band width 
will only be noted briefly In order to exemplify pertur- 
bations to the case where Eq. |T]) only includes a discon- 
tinuity, results will often be visualized in the case where 
only the first two terms in the sum of Eq. ((T|) are present, 
i.e. p(e) = (p + pie) ■ 6(e). 

From this density of states the doping x can be analyt- 
ically calculated for a given chemical potential p as 



de(l -f(e-p))p(e) 



(2) 



where fix) represents the Fermi function. The integra- 
tion occurring in Eq. @ in combination with Eq. ([1]) can 
be performed using polylogarithm functions^: 



$>/r (ri+1) 

n 
n 

3=0 



1 

7 Lij+i ( 

in ~ 3)1 



u{fi)e 



-pw\ 



dt 



(3) 



where the polylogarithm functions 
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n\ J xexp(r) — 1 

as well as the renormalised values and the negative fu- 
gacity, 
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were introduced. The exact solution of Eq. ^ would give 
the temperature dependence of the chemical potential. 
Unfortunately, it can only be performed numerically, but 
useful insight can be gained by different approximation 
schemes. Those will be discussed in the following section. 

III. THEORETICAL BACKGROUND 



III. A. Sommerfeld Region 

From Equation §5§ the result of Sommerfeld expansion 
can easily be reproduced: By expanding the polyloga- 
rithm functions for large negative arguments^, respec- 
tively low temperatures (/3p 3> 1), up to the first polyno- 
mial terms in temperature, 

L"/ 2 J (_o..\n-2m 

Li„(-e"^) = 2 LW-1) \ m 9 „ 
^— ' [n — 2m)! 

m=0 y ' 



(-l)"Li„(-e^) 



(5) 



6 n\ 

0{{Pp) n - A )+0{e^), 



(6) 



a quadratic equation for the chemical potential is ob- 
tained, if the terms of order 0(/3~ 4 ), ©(e^) and the 
effects of a finite band width W are neglected. Expand- 
ing its solution for small temperatures yields the familiar 
result of Sommerfeld expansion: 



k B l T F w ■ [ — 



(7) 



where the Fermi liquid parameter vj and the Fermi tem- 
perature Tp are given by 



7T 2 2 p'(T F /h 
— "TT^B 

6 



p(T F /k B ) 



and Tp 



y/1 + 2yn - 1 
k B ri 



(8) 

with p' (e) denoting the first derivative of the density of 
states. For the latter result was assumed that the coef- 
ficients in the Taylor series, Eq. ([1]), are decreasing with 
growing order, i.e. r n <C T\ for n > 1. 

When Sommerfeld expansion is applied to determine 
the low temperature behaviour, the convergence of the 
Taylor series used therein is usually not considered. For 
a Fermi energy close to a discontinuity this induces a 
further breakdown criterion. It can be obtained when 
recalling that the expansion (JSJ is made in the case (3p 3> 
1. As a measure of the condition where this expansion is 
violated, the temperature where (3p = 2 can therefore be 
taken. With this the criterion the expansion ^ shows 
that the Fermi liquid expression should be valid up to the 
temperature 



y/1 + Tpw/k B - 1 

w/k B 



(9) 



Therefore the neglected terms in this approximation ap- 
plied to Eq. (U]) are smaller than v~ x = — cT 2 . While the 
temperature Tg will always be lower than the Fermi tem- 
perature Tp, the comparison to the numerical solution of 
Eq. in Fig. [2^ shows that it seems to be important 
only when the effect of the discontinuity in the density of 
states is more important than that of the slope. Accord- 
ingly, in the case of a dominating discontinuity Ts de- 
scribes the numerically studied breakdown of the Fermi 
liquid more accurately. In addition, the breakdown tem- 
perature T$ and the Fermi temperature Tp as well as the 
numerical solution show for large renormalised doping y 
and small parameter ratios r n a region where Sommer- 
feld expansion is valid at room temperature, and even 
beyond. 

III.B. Boltzmann Region 

On the opposite side of physical temperatures, the re- 
sult of Boltzmann approximation can be deduced as a 
Taylor expansion of the polylogarithm functions in v in 
the limit (3p <C — 1. The result in first order is given by 



Ms 



-k B T In 



k B Tf{k B T) 



(10) 



where the quantity r(e) = ~^2 n r n e n was introduced 
which will be discussed in the following section IIII.CI 
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Figure 2: Relative errors of the chemical potential obtained by the discussed approximations to the one calculated 
numerically from Eq. as functions of the parameter ratio r± and the renormalised doping y for a specific temper- 
ature of T = 400 K and a density of states given by p{e) = (po + p\e) 6(e). The shading in Fig. [2^ runs from 0% 
(darkest) to 10% (lightest blue) for the comparison of the solution from Sommerfeld expansion Eq. ([7]), and from 0% 
(darkest) to 10% (lightest orange) for the Boltzmann approximation Eq. (|10p . respectively 0% to 100% for the one 
using the APLD result Eq. (|2ip in the inset of Fig. [2b. The main figurc[2b shows the relative errors for the quadratic 
expansion of the APLD Eqs. (|3"3")) and with the used expansion temperature T — 800 K. Additionally drawn 
are the curves where the specified temperature T matches the Sommerfeld breakdown temperature Ts, the Fermi 
temperature Tp, the Boltzmann breakdown temperature Tb, the degeneracy temperature Td og , the temperature 
Tp, where the APLD has a logarithmic singularity, or the transition temperature To, where the chemical potential 
equates with the negative thermal energy:/! = — k B Tg. In the inset of (b) the degeneracy temperature Td eg as well 
as Tp are very close to the points where the relative error between the APLD and the numerical solution increases 
rapidly. 



In case the band width W is taken into account the 
defined quantity is shifted in this approximation to 
f(e) — r n £ n T(n + 1, W/s)/n\, where T(n,x) denotes 
the upper incomplete gamma function. Since this func- 
tion becomes an exponentially decreasing function for 
large second argument, i.e. T(n,x) — > e~ x for x — > oo, 
the effects of a finite band width can be neglected if the 
thermal energy k B T is still smaller than the band width. 

In case the band width is neglected a similar break 
down temperature as in Eq. @ can be extracted for 
Eq. ffTTJf) but for the criterion f3p = —2. This means the 
expression (|10[) should be valid for temperatures larger 
than 



Tp, 



i/l + 4yr ie 2 
2fc B ri 



(11) 



where corrections from this approximation applied to 



Eq. ^ are of order v 2 



4 and where again a decreas- 



ing importance of the Taylor coefficients with growing 
order was assumed. 



will then show a finite linear term in contrast to Eq. |7|): 



Mb = /4 0) + /4 1} ■ (T — T) + -pf ■ (T - f f with 
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k B fr(k B f) 
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(12) 
(13) 



_fc Bln MW) -k B [l + k B f (lnf)«], (14) 



-^[1 + 2k B f (Inr) (1) + (fc B T) 2 (lnf)( 2 )] 



(15) 



where (\nf)( n > denotes the n-th derivative of the loga- 
rithm of the quantity f at the expansion energy k B T. 
These terms are not described within the Fermi liquid pic- 
ture, establishing instead a region of an apparent Fermi 
liquid (AFL). In particular, for an AFL behaviour a non- 
vanishing linear term in the chemical potential is ob- 
tained which even does not vanish when reordering the 
coefficients like in the Fermi liquid case 
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(/4 0) 



pff 2 /2) 
{p^-2p%T)-T + p%/2-T\ 



In order to bring the expression (|10p into a similar form 
as obtained in the low temperature expansion it can be 
linearised around a certain temperature T . However, it 



However, due to the expansion around a certain tem- 
perature the coefficients depend on this quantity, too. 
Nevertheless, for large temperatures the terms in squared 
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brackets in Eqs.[*Hl and [Tulare constant in leading order, 
resembling only a weak dependence of the coefficients. 

III.C. Intermediate Region 

In order to fill the gap between Boltzmann and Som- 
mcrfcld region, it is useful to take a closer look at the 
similarity of the poly logarithm of different orders. To 
that aim it is useful to recall that every polylogarithm be- 
haves similar at small negative arguments, e.g. Li(0) = 0, 
Li'(0) = 1. Since the first polylogarithm resembles an 
elementary function Lii(z^) = — ln(l — v), it is very ap- 
pealing to consider the other orders in terms of an ap- 
proximation of the polylogarithm difference (APLD) 



d n {v) = Li„ + i(i/) - Lii(^) . 



(16) 



occurring in Eq. ^ , in terms of a Taylor series around a 
negative expansion fugacity vq: 



(I 



/r V~ (/r 1 ) in(i - v) - r n <Uy ) /r (n+1) 



(17) 



where for simplicity again the limit W — > oo was consid- 
ered. 

The first terms, proportional to the logarithm, in this 
approximation can be understood as following from a 
temperature dependent, averaged density of states. Re- 
viewing the expressions of the previous subsection it has 
already occurred as the denominator of the fugacity in 
the result of Boltzmann approximation in Eq. (|10p . In 
the case of a dominant discontinuity or slope it takes the 
form 

r( S )np(k B T)/p , (18) 

but provides different expressions for more general den- 
sity of states. Additionally, the zeroth order of the Taylor 
series only shifts the doping value used in the theory 



y -> y + dy(k B T) with dy(e) = 2J r n d n (v ) £ 

n=£0 



n+1 



(19) 

since this term does not depend on the chemical potential. 
In contrast, the first order terms in this expansion would 
result in the chemical potential given by 

k B Tr{\ + e^) exp(fc B Tr (1 + e^)) 

u rr- ( V + dy{k B T) - k B T r (up -1) \ 

= fcBTrexp l wwo J' (20) 

where k B Tf is the prefactor of the linear term in Eq. Q17p 
similar to dy{k B T) for the zeroth order. The solution 
of this equation can be found as the Lambert W func- 
tion— with the right-hand side as argument. Of course, 
for small (positive) arguments this function can again be 
linearised. Since this leads to the same solution as if only 
the zeroth order term is taken, the first order contribu- 
tions only give corrections at large doping, Eq. (TT9"j) . or 



low temperatures. Therefore this as well as all higher 
orders will be neglected in the following. The chemical 
potential in this approximation is then given by 



fip = k B T In 



exp 



y + dy{k B T) 



k B Tr(k B T) 
k B T (z + \n(2sinhz)) , 



- 1 



(21) 
(22) 



with the argument of the hyperbolic function 

_ y + dy(k B T) _y + E„^o r n dM (k B T) n+1 



2k B Tf{k B T) 



2 En'"" (fc B T)«+l 



(23) 

The validity of this expression can be related to the 
results of Boltzmann approximation and Sommerfcld ex- 
pansion in particular cases: For a very large expansion 
temperature the difference of the polylogarithm d n (fo) 
becomes negligible. Due to the assumed high temper- 
atures the result of Boltzmann approximation \ip in 
Eq. pop is obtained by an additional expansion of the ex- 
ponential function in Eq. (|2ip . Contrastively, in decreas- 
ing the expansion temperature the fugacity will reach 
unity, i.e. Vq — > — 1. After this point a slight variation 
of the chemical potential or the temperature will lead to 
a greater change in the fugacity. Therefore higher or- 
ders become important and the Taylor scries in Eq. (JT7J) 
will break down for smaller temperature differences. This 
means the APLD is only appropriate for temperatures 
larger than the degeneracy temperature Tdcg where the 
chemical potential vanishes. According to Eq. ([H]) this 
temperature is implicitly given by 



y 



E r ™ Li„+i(-l) (k B T deg ) 



n+1 



(24) 



If the Taylor series in Eq. (PTJ) is dominated by its first 
terms, it can be solved yielding 



T, 



deg 



In 2+ v/(ln2) 2 + ttVi/3 
Tr 2 k B ri/6 



(25) 



However, for a very dominant discontinuity, the first term 
in Eq. (|17p will give the dominant contribution to the 
doping. Therefore, the other terms will resemble small 
corrections whose importance will decrease as the first 
term increases. Since the latter was treated exactly, this 
means that the validity region is increased to very low 
temperatures. In particular, for the case the density of 
states is only described by a discontinuity the APLD be- 
comes exact and resembles the result of Sommcrfeld ex- 
pansion after performing a Taylor series for small tem- 
peratures. 

Surprisingly, at finite values of the ratios r n and the ex- 
pansion parameters d n {vo) the chemical potential in this 
approximation has a logarithmic singularity where the ar- 
gument of the exponential function in Eq. (|2ip vanishes. 
In order to avoid this divergence Boltzmann approxima- 
tion should be used for temperatures similar or larger 
than this temperature Tp which is given for the first two 
terms of the Taylor series in Eq. {1} by 



y 



-di(v )rik B < 



(26) 
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A suitable expansion point which lies between the Som- 
merfeld and Boltzmann breakdown temperatures given 
by Eqs. @ and pip is the transition temperature To 
where the exact chemical potential, as given by Eq. ([3]), 
is equal to the negative thermal energy p = —k B T or 
vq = — e _1 . Similar to Eq. (|25p this temperature can be 
stated explicitly if the Taylor series in Eq. fT]) is domi- 
nated by the first terms. In this case it is given by 



To 



ln(l - v Q ) - \J (ln(l - v )) 2 - 4yri Li 2 (v ) 
2k B n Li 2 (^ ) 



(27) 



With this expansion temperature To further simplifica- 
tions can be brought to Eq. (|2ip : In this case the absolute 
value of the difference of the polylogarithm functions of 
first orders amounts to |di(fo)| = 0.03 and monotonically 
decreases with larger expansion temperatures. Multiply- 
ing this number by the typical values for k B 2 T 2 at room 
temperature leads to a result which is much smaller than 
the renormaliscd doping y. Although the difference in- 
creases slightly for higher orders, the shift of the doping 
dy in the result of the APLD, Eq. (|2"Tj) . may be neglected, 
especially if the Taylor series is governed by its first terms. 
Since this means neglecting even the first order in the ex- 
pansion of the APLD (|17|) the expansion point i^o would 
not be present in the solution of the chemical potential 
Eq. (|2ip . Therefore an even larger application region is 
expected in this case which is in agreement with the di- 
verging behaviour of the break down temperature Tp as 
well as with the discussion of the lower break down tem- 
perature after Eq. (|2"5j) . 

Since in this limit the results of the APLD can be 
understood as those following from a constant averaged 
density of states given in Eq. (fT5|) . at first the chemi- 
cal potential of only a discontinuity will be studied, i.e. 
p(e) = po <d(e)@(W - e). In this limit, where the APLD 
becomes exact, a finite bandwidth W shifts the chemical 
potential by 



/ip — > pp — k B T In 



1 — exp 



V 



w 



k B Tf(k B T) 



(28) 



which resembles a similar contribution as Eq. [2T] but for 
the upper band edge. Following from this shift the chem- 
ical potential reaches a constant value in the high tem- 
perature limit j3 — > 0, 



p P 



k B T In 



V 



W-y 



(29) 



in contrast to the chemical potential obtained by Boltz- 
mann approximation in Eq. (|10p . Therefore, the assump- 
tion of that approximation /3/x <C — 1 is not valid when 
the thermal energy surpasses the band width W . How- 
ever, for a sufficiently large band width W 3> k B T + y 
with respect to the temperature one is interested in, the 
shift can still be neglected. 

In order to compare the observed behaviour to the one 
known from a Fermi liquid as derived in Eq. Q, the 
chemical potential of Eq. (|2"2"]) can be brought into a sim- 
ilar form. Performing a quadratic expansion around the 




Figure 3: Dependence of the coefficients of the quad- 
ratic expansion of the chemical potential in APLD 
/Up = ip — q e SoT + zuT 2 on the expansion temperature 
T for fixed renormalised doping y = 0.2 eV and parameter 
ratio n of different regimes. For small parameter ratio 
ri as given by Eq. (|30p (thick lines), the temperature 
dependence of these parameters are weakest around the 
temperature Tp ^ given by Eq. (|3ip (dashed vertical line). 
For significant parameter ratio T\ = 2 (eV) -1 (r n = for 
n > 1) according to Eqs. (|33p and (|34p (thin lines), the 
stationary points have increased and have begun to shift 
independently, causing the AFL to destabilize. In the 
figure ip, q e So, w are given in units of fee -1 and 
yk B ~ 2 respectively 



temperature T leads to 



& = pf(f) + • (T - T) + -pf(f) ■ (T - ff 



= -k B [z + ln(2 sinh z)] T 

+ k B [z coth z - ln(2 sinh z)\ (T - f) 
z 2 1 



' 2 sinh 2 z f 



- (t - ry 



(30) 



with z — y (2/cbT) -1 , the argument of the hyperbolic 
functions as in Eq. (|2l?l) . Due to the fact that the linear 
coefficient after reordering will not vanish, which can be 
seen in Fig. O this leads to an AFL behaviour. Further- 
more, since the quadratic coefficient ^pp(T) possesses 
in its T dependence an inflexion point near its maximum 
at 

3 

zcothz = - , (31) 

the terms of higher order will have small values near this 
maximum. Therefore the quadratic expansion will ap- 
proximate the function very well around this point defin- 
ing a particular expansion temperature. Since the numer- 
ical solution of Eq. (j3"Tj) reads z f» 1.29 this temperature 
is given by 

)5— -y = 0.8 T s = 0.3T dog . (32) 



1 p 



eV 



An interesting aspect of this temperature comes from 
the fact that the chemical potential have not yet moved 
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far from the Fermi energy. In fact, the low temperature 
result obtained in Eq. (J7|) would describe a chemical po- 
tential independent of temperature. It therefore remains 
at the Fermi energy ef with increasing temperature. At 
the specified temperature, when it differs from the Fermi 
energy by (ef ~ (J,p(T p ')) /ef — 3% it starts moving 
towards the band edge. This is due to the fact that in 
the doping constraint, Eq. the Fermi distribution is 
cut at the band edge where it is reduced by 8% at T p '. 
This start of the motion of the chemical potential can be 
approximated as the Taylor series in Eq. (|30p with non- 
vanishing quadratic coefficient, therefore distinguishing it 
from the Fermi liquid parameter. Nevertheless, a cross- 
over from a true Fermi liquid behaviour, i.e. with finite 
parameters, to such an AFL can not be observed in this 
limit. 

However, on the one hand in considering very small 
additional coefficients of higher order terms in the Tay- 
lor series, the real Fermi liquid parameter would be non- 
vanishing. On the other hand, the boost of the chemical 
potential near the temperature T p may still be possible 
to observe on top of the background behaviour caused by 
these coefficients. Since the breakdown temperature of 
the observable Fermi liquid T$ will be slightly lower than 
the temperature where the AFL behaviour occurs, this 
leads to a comparably small transition region between a 
Fermi liquid and an AFL behaviour. Analytically, this 
represents itself as follows: When taking the effects in 
Eq. (|22|) of the full Taylor-series of the density of states 
into account, both the hyperbolic argument as described 
in Eq. (|23p and the modification of the derivatives in the 
quadratic expansion have to be taken into account: 

Hp\f) -> n ( p\f) - k B z (1 + coth z) , (33) 

/ z\ 2 l + cothz 

zd — > zd 1 H — + k B ^ 

V z) 2(f(k B T)) 2 

■ [zk B ff {2) r(k B f)+dy {1) f {1) 

-2zk B f (rW) 2 - l -dy^ r(k B f)} , (34) 

with z = (2zk B ff^ - dyM)/(2f) and where ft") de- 
notes the n-th derivative of the renormalised density of 
states and dy^ respectively the one of the doping change 
both taken at the thermal energy k B T. These adjust- 
ments cause the common stationary point of all three 
coefficients of the chemical potential /ip', as a polyno- 
mial in T, to shift independently. As can be seen in 
Fig. [5b, this destabilizes the quadratic approximation 
thereby destroying the AFL, stressing the role of a large 
discontinuity for this kind of AFL. However, in adjust- 
ing the renormalised doping y to given parameter ratios 
r„, which should not be too large, a narrow region re- 
mains where the different stationary points are close to 
one another. As indicated in Fig. [3] the observation tem- 
perature Tp 9 "* of the AFL of this region would then be 
larger than the temperature given in Eq. (|32p . Since the 
degeneracy temperature given by Eq. decreases with 
growing parameter ratios r n , this even leads to an AFL 




<1 [(eV)" 1 ] 

Figure 4: Phase diagram for T = 400 K (areas) and 
T = 800 K (lines). In the different coloured regions the 
relative errors between the exact solution of the chem- 
ical potential and the one described by the stated be- 
haviour is smaller than 50%. For the AFL by APLD 
(AFL region on the left), Equations ([3"3")) and (|34j) arc 
compared to Eq. (|21j) . respectively for the one given by 
Boltzmann approximation (AFL region in the lower right 
corner) Eq. (|I2p to Eq. pup . In these regions the more 
general Fermi liquid and AFL expressions ([7]), (fTT)|) , (j2Tj) 
differ from the exact solution by less than 1%. 

behaviour where the chemical potential has left the band 
as can be seen in Fig. [U 

Reviewing the results of this section characteristic 
(breakdown) temperatures were found in Eqs. ([H]), ©, 
P~Tj) . (|2"5 ]) -([2T P which may be ordered as 

T s < Tdcg < T < T B < T P . (35) 

While for large discontinuities Tp lies between the Som- 
mcrfcld breakdown temperature Ts and the degeneracy 
temperature Ideg it becomes larger than T^ eg for small 
ones. In contrast to the AFL by Boltzmann approxima- 
tion, where the temperature has only to be larger then 
the Boltzmann breakdown temperature Tp, the AFL 
by APLD is centred around the temperature T p ' from 
Eq. (|32p for vanishing parameter ratios r n or slightly 
above for finite values of it, and crucially depend on the 
influence of the discontinuity. This can be seen in the 
(ri,j/)-phasc diagram as found in Fig. 3J Of course in a 
more realistic density of states the presence of an upper 
band can destroy these phases if the chemical potential 
becomes too close to itQ. Therefore a large enough gap is 
necessary, too. 

IV. APPLICATION 



corresponding Taylor series of an upper band separated by a 
gap A would lead in Eq. Q similar to a lower bandwidth to poly- 
logarithm functions with argument ue^^. Therefore the terms of 
the approximations would be corrected only by exponentially sup- 
pressed terms which are negligible if the chemical potential remains 
smaller than the size of the gap. 
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In this section the formulas obtained above are applied 
to doped CuCr02 and doped CuRh02- Since obtain- 
ing the chemical potential experimentally is a challenging 
task, the temperature dependence of the thermopower 
should be used to determine this potential within the 
TICR framework and, with the results of the last sec- 
tion, an effective density of states. However, fitting an 
infinite number of parameters as appearing in its Tay- 
lor series, Eq. ([T]), is obviously not practicable. Instead, 
only the first two terms of the Taylor series should be 
considered. In addition, the replacement of chromium by 
magnesium as well as the substitution of rhodium later 
leads to a hole-doping in the t^g band. For this kind of 
doping the upper band edge is important in contrast to 
the lower one which was described in Eq. ([IJ for electron- 
like doping. In order to apply the previously discussed 
framework, the signs of all energies should therefore be 
inverted. This leads to the density of states to be given 
by 

P (e) = ( Po - Pl e) ■ B(-e) . (36) 

Furthermore, it will turn out, that the effective charge 
carrier density is accessible, too. 

IV.A. The TICR framework 

The thermopower can be studied within the approach 
of the temperature independent correlation functions ra- 
tio approximation (TICR) 12 i 25 . This approximation is 
derived from the Kubo formula where the thermopower 
is given b y 18 ' 43 



S = 



1 



(jEjn) 
<leT V (jnjn) 



M 



(37) 



where q e is the negative charge of the electron and the 
particle (energy) current operator is denoted by j n (js)- 
If a shift of the single-particle energy is considered weight 
can be shifted between both terms on the right side of this 
formula. In particular a shift Eq exists, where the trans- 
formed correlation function between j n and je vanishes 



(j' E j'n) = (jBjn) + Eo(jnjn) = . 



(38) 



Therefore the thermopower is only given by the trans- 
formed chemical potential fx' = fx — Eq: 



S=-^{E a -fx) 

ye 1 



(39) 



Comparing this result to Eq. (|3"T)) determines the shift 

O'sin) 



En = 



\3n3ri 



(40) 



which would be temperature dependent in general. Thus 
although Eq. (|39|) describes the transport function S in 
terms of the equilibrium property fx' , for the full temper- 
ature dependence transport theory is still necessary to 
determine the shift Eq. 

If the transport through a single level is considered, 
this shift can be determined in case the coupling to bal- 
listic leads is weak and temperature independent. If such 



a level is in resonance, i.e. at the same energy as the 
chemical potentials of both leads, a temperature gradi- 
ent would not lead to a current / due to the symmetry 
of the model. This means that the thermopower would 
vanish. Thus in the following will be considered the oppo- 
site situation where the level would be far away from the 
chemical potentials. In this limit, the Lorentzian trans- 
mission T{e) for the non-resonant level model might be 
expanded in terms of the broadening Y of the level 



T{e)=T8{e~eQ) + 0{T 2 ). 



(41) 



where Eq denotes the energy of the level and the normal- 
isation at zero temperature was taken into account. The 
Landauer formula 44 ' 45 for the tunneling current / gives 
for such a transmission in the spinlcss case 

1 = S / r(e) " ml) " Ir{£ - de (42) 



irk 



[/l(£o - Ml) - Ir(so - Mfl)] 



(43) 



where H is the reduced Planck constant and the Fermi 
function / as well as the chemical potential fi in the left 
and right lead where distinguished by the indices L and 
R as will be later on for the temperature as well. For the 
condition that no current flows / = 0, as assumed in the 
definition of the thermopower, the Fermi functions and 
therefore their arguments have to be equal 



(eo - Mr) ^bTl = k B T R (e - fx L ) 



(44) 



In expressing therein the quantities of the left and right 
environment by their small differences q e AV and AT as 
well as the averaged temperature T and chemical poten- 
tial fx, the thermopower is obtained as 



S = 



AV 



1 



AT q e T 



£o 



M) 



(45) 



Therefore the Eq function assumes the constant position 
of the considered level £o. Furthermore, the same rela- 
tion can be derived when the transmission is governed by 
thermal processes^. In addition, it can be shown that 
this shift is not important at very high temperatures 1 ^. 

Motivated by these results the basic assumption of 
TICR treats the shift Eq independent of temperature 
in the considered temperature region. Especially when 
the chemical potential is inside the band gap where the 
structure of the bands should not be important, the sit- 
uation is similar to the case of the non-resonant level 
model. Therefore the approximation should be valid in 
this case. The TICR can thus be viewed as a general- 
ization of a purely stochastic approac h 18 i 39 i 40 where the 
first order corrections in the broadening T is taken into 
account. However, if the thermopower for a usual three- 
dimensional dispersion relation is calculated in Fermi liq- 
uid theory—, it deviates from the result given by TICR 
for the chemical potential following from such a disper- 
sion relation. This means the TICR will break down at 
low temperatures for a metal when the chemical potential 
is inside the band. 
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In contrast to earlier studies where the shift Eq was set 
to the Fermi energ y 35 i 38 or omitte d 18 i 39 i 40 it has proven 
useful in recent studie s 12 i 25 i 26 to keep it as a free param- 
eter, especially at large temperatures. Furthermore, this 
constant may account for some of the contributions due 
to scattering processes caused by disorder which might 
be important at these temperatures. Additionally, it pre- 
serves the freedom of choosing an arbitrary zero-point of 
the single-particle energies. In the following investigation 
the stated values of the constant will be with respect to 
the choice made in Eq. ([56")) while the chemical potential 
is given by the expressions determined in the last section. 
If, as found there, a chemical potential of an apparent 
Fermi liquid (AFL) is assumed, 



/i = ef — q e SoT + wT 



(46) 



the thermopower will be T-linear with a hyperbolic and 
a constant offset: 

S{T) = ^-^ + Sa + ^-T. (47) 

Qe 1 \Qe\ 

It therefore can describe the observed linear thermopower 
with an extrapolated offset So as shown in Fig. Q] if the 
hyperbolic behaviour can be neglected due to high tem- 
peratures. 

From the expression (|4"7)l it is easily seen how to deter- 
mine the linear and quadratic coefficient So, w of the 
chemical potential out of the temperature dependence 
of the thermopower. However, the TICR parameter Eq 
makes it difficult to determine the constant offset. While 
at very large temperatures one can try to neglect the Eq 
constant, 

|£o|«|e>|, (48) 
at finite temperature the study of the quantity 

TS{T)-fS{f) 
T-f ' 



(49) 



where T denotes a reference temperature, might be use- 
ful since it does not depend on the TICR constant Eq 
and should be an exact linear function with respect to 
temperature. The experimentally accessible parameters 
in Eq. (|47p allows to obtain the parameters of the effec- 
tive density of states introduced in Eq. (|36|) and renor- 
malised in Eq. ^j. However, their values will depend on 
the different approximation schemes discussed in Sect. Hill 
and therefore labeled below by an index. Furthermore, 
assuming common parameters of the effective density of 
states for all samples of an experimental series, e.g. for 
several manganites^, the exact effective charge carrier 
density can be obtained in relation to a reference sam- 
ple. Therefore this technique resembles an alternative 
treatment to identify the doping values of such a series 
complementary to other experimental based techniques 
like EDS. 

Combining the TICR expression (|59"|) with the result of 
the APLD given in Eq. ([2"9"]) for only a discontinuity in the 
density of states in the limit of large temperatures yields 
the famous Hcikes formula in the atomic non-degenerate 



limit poW = 1, or modified version of it like derived ini£ 
with g = poW as the degenerate factor. Furthermore, the 
combination of the TICR with the result of Boltzmann 
approximation in Eq. (|10j) leads to the observation, that 
the thermopower still resembles the doping dependence 
of this formula at low doping S ~ — lnx. The prefactor 
can also be related to the findings reported for the ther- 
mopower in the atomic limit of degenerate energy levels. 
For a lattice consisting of atoms with g fully occupied de- 
generate states a doping with atoms missing one fermion 
would lead to a g-fold degeneracy of the ground-state of 
the latter atoms. As discussed in— this would lead to a 
thermopower of 



a - kB i ( 1 ~ x 

•^atomic — m 



(50) 



in the high temperature limit. When fcrmions dclocalisc 
the distribution in the density of states broadens. How- 
ever, since the total number of states remains fixed a 
discontinuity in a density of states like Eq. ([56]) of 



Pa = 777 



w 



(51) 



is necessary if the slope therein is not considered. Ac- 
cording to Eqs. p2[) and (|47[) the thermopower will be 
given by 



S=^ln 



k B T 1 



q e 2T 



T. 



(52) 



if the hyperbolic terms is omitted due to high temper- 
atures. While for large expansion temperatures T the 
linear term will be negligible the degeneracy factor of the 
constant term is renormalised compared to the atomic 
limit. Therefore in considering a broad band only the 
thermally accessible region of the band enters in the ther- 
mopower. 

IV.B. CuCn x M gx o 2 

Carrying out the analysis for the strongly correlated 
delafossite CuCri_ x Mg x 02 gives already reasonable pa- 
rameters of the density of states using Boltzmann approx- 
imation for large parameter ratio, in contrast to Sommer- 
feld expansion. Furthermore, the results coincide with 
the one obtained by the APLD if a larger expansion tem- 
perature To — > p{2k^)~ 1 is used indicating a very low 
temperature scale. This requires to take the first order 
in this approximation proportional to dx(vo) in Eq. (JTTJ) 
into account, too. 

While the results of the TICR parameter Eq = 
100 meV are close to the ones reported for doped man- 
ganites, but with opposite signs since they are electron- 
doped—^, the slope p\ seems to dominate the effective 
density of states for this material. Although this is in 
agreement with a previous first principle study—, it chal- 
lenges the determination of the small discontinuity pQ. 
Furthermore, the results exhibit a strong increase of both 
parameters for the highest doped samples which falls in 
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Figure 5: Comparison of the thermopower of 

CuCri_ x Mg x 02 given by experimental data for nominal 
doping x = 0.5% (stars), x = 1% (crosses) and x = 2% 
(circles) to the one obtained from the parameters of the 
density of states Eq. (|53"j) for values of the effective charge 
carrier density of 0.35% (solid line, with Eq = 134 meV), 
0.85% (dashed line, with Eq = 137 mcV), and 1.1% (dot- 
ted line, with Eq = 105 meV). Inset: Temperature de- 
pendence of the chemical potentials for the corresponding 
parameter values. A quadratic behaviour can be observed 
at very low and very high temperatures. 

the range where the formation of a secondary spinel or 
CuO phase was previously reporte d 11 ' 48 . They are there- 
fore omitted in the following investigation. The other 
samples show still results of the parameters of the effec- 
tive density of states depending on the doping in contrast 
to the parameter ratio n = 200 (eV)" 1 . This encourages 
adjusting the doping to find the effective charge carrier 
density, as was done for several other manganites— . For 
this material this problem was experimentally addressed 
in a previous publication™ leading to quite different dop- 
ing values by statistical electron diffraction spectroscopy 
(EDS). If the most reliable value of the reported data 
x = 1.1% is used for the nominal one of x = 2% the 
parameters of the density of states, 

Po.b = 0.121 (eV)" 1 and p lmB = -24 (eV)" 2 , 

(53) 

coincide for all considered samples for effective charge 
carrier density of 

x = 0.35%, 0.85%, 1.1%, (54) 

which arc within the error margins of the one suggested 
by the EDS measurements. Due to this renormalisation 
the Eq parameter are only slightly modified during the 
minimization process. 

The calculation of the thermopower using the exact nu- 
merical solution leads to good agreement between theory 
and experiment as can be seen in Fig.0 if the TICR con- 
stant Eq is adjusted to the modified doping values by an 
additional fit. Additionally, the hyperbolic offset caused 
by the TICR parameter seem to describe the increase at 





x — 1% 

,, „( _ ) 
fip fl p fl p 


x — 4% 
fip Hp fip 


x = 10% 

,,(?) ,,(-) 
fip fip fip 


ri [(eV)- 1 ] 
y [mcV] 
E [mcV] 


3.53.1 - 
4.84.6 2.4 
33 31 6.1 


3.4 3.4 3.4 
6.0 6.0 3.1 
11 11 -14 


3.4 2.5 - 
21 19 13 
36 32 17 



Table I: Stationary points of the fitting routine when 
adjusting the parameters in the formulas of the APLD in 
a region from 600 K to 1000 K to the experimental data 
of CuRh 1 _ x Mg x 2 . 



low temperature and doping, while the agreement is lost 
below T ~ 600 K for intermediate doping values. This 
might be due to a breakdown of TICR, but surprisingly 
coincides with the Boltzmann breakdown temperatures 
for these samples T B = 356 K, 497 K and 608 K. These 
values confirm the placement of the observed behaviour 
in the AFL by Boltzmann approximation, too. 

IV.C. CuRIn x M gx 2 

The application of the presented technique to the mea- 
surements of the doped band insulator C11RI1O2— using 
established approximations does not give cohesive results. 
In comparing the results stated in Tab.|T]and produced by 
a fit using the formula of the APLD with the expansion 
temperature given in Eq. (|27p . the values of the param- 
eter of the TICR Eq ~ 30meV are in agreement for the 
lowest and highest doped samples and fall in a region 
previously reporte d 12 i 25 , too. Furthermore, the results 
from the expansion for small parameter ratio obtained 
in Eq. (|30p give values similar to the ones obtained from 
the quadratic expansion and from fitting the full solution 
stated in Eq. (|2"Tj) . Again the found parameter ratios are 
close, while the values of the effective doping vary. There- 
fore the doping can be adjusted where the results of the 
parameters of the density of state coincide for values of 

x = 2.3%, 2.8%, 10% . (55) 

With these values of the doping all samples give therefore 
rise to the same parameters of the density of states of 

p ,p = 4.7 (eVr 1 and p hP = -16 (eV)^ 2 . (56) 

As can be seen in Fig. [(JJ the results of the APLD are 
close to the ones obtained by first principle studies. The 
extracted parameters places the largest doped sample of 
this material in the phase diagram shown in Fig. Q] pre- 
cisely in the narrow area (intermediate parameter ratio 
n) of the validity range of the AFL by APLD, while the 
lowest doped samples are closer to the region described 
by Boltzmann approximation. This is further supported 
by the fact that the Boltzmann breakdown temperature 
for the largest doped sample Tg(x = 10%) = 1316 K is 
larger than the one considered while the expansion tem- 
perature of the APLD T (x = 10%) = 653 K lies within 
the region where the linear behaviour is observed. There- 
fore the observed AFL can be understood as the remains 
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Figure 6: Density of states of CuRhi_ x Mg x 02 as given 
by the augmented spherical wave 2 method (crosses) and 
their shape as determined by a fit to the exact form 
Eq. (|21[) (thick) and quadratic expansion Eqs. Q33 |) -(|34 p 
(thin) in APLD for different values of the doping of 
x = 1% (solid), 4% (dashed) and 10% (dotted). 



of a larger manifestation for larger discontinuity. Com- 
pared to doped CuCr02 this material exhibits therefore 
a discontinuity enhanced thermopowcr although larger 
doping values were considered. Furthermore, the recal- 
culation of the thermopowcr using the exact solution 
shows extraordinary good agreement as can be seen in 
Fig. [7] At low temperature the hyperbolic offset Eq of 
the thermopower dominates the theoretical behaviour for 
this material, too. Nevertheless, since the data does not 
show a strong upturn this rather indicate the breakdown 
of the TICR. However, a slight upturn can be found for 
the lowest doped sample which compared to the inter- 
mediate one has a larger Eq constant. Since interaction 
effects are expected to be negligible due to nearly filled 
bands this might arise from scattering processes caused 
by a larger disorder in this sample which could explain 
the insulating behaviour at very low temperatures^, too. 

From the parameters of the density of state further 
quantities like the specific heat can be determined which 
is given in Fig. [SJ Since the determination of the specific 
heat is not dependent on the TICR parameter a Fermi liq- 
uid behaviour can be observed at very low temperature. 
Nevertheless, for low doping the effective mass gets renor- 
malised in comparison to calculations in which the discon- 
tinuity is neglected. This phase breaks down at a temper- 
ature which decreases with decreasing doping, in agree- 
ment with the expected insulating behaviour for zero dop- 
ing. However, the breakdown temperatures are obviously 
not well described by the Fermi temperature Tp, but 
seem to occur at temperatures somewhat below the ana- 
lytically introduced smaller breakdown temperature Tg. 
In the vicinity of this temperature a peak in the deriva- 
tive of the specific heat is found, as shown in the inset of 
Fig. El It should be observable in experiments, in partic- 
ular for low doping. Above the transition area, around 
the degeneracy temperature T^ es , a linear region emerges 
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Figure 7: Comparison of the thermopower of the de- 
lafossite CuRhi_ x Mg x 02 given by experimental data for 
nominal doping x = 1% (stars), x = 4% (crosses) and 
x = 10% (circles) to the one obtained from the param- 
eters of the density of states Eq. ([55]) for values of the 
effective charge carrier density of 2% (solid line, with 
E Q = 30mcV), 3% (dashed line, with E = 15meV) 
and 10% (dotted line, with E = 36meV). Inset: Tem- 
perature dependence of the chemical potentials for the 
corresponding parameter values. Two quadratic regions 
can be identified, one described by Fermi liquid theory 
at low temperature and one in the framework of an AFL 
at higher temperatures similar to the thermopower and 
conductance shown in Fig. [TJ 



which extends to Tq, where the chemical potential has 
left the band, and for large doping even beyond. Note 
that the difference between the apparent effective mass 
and the bare one exceeds the one usually expected from 
interaction effects obtained using the slave-boson saddle- 
point approximation applied to a nearly filled band—. At 
higher temperatures for the two lower doped samples, the 
slope gradually reaches a slightly smaller value, forming 
another linear region above the breakdown temperature 
of Boltzmann approximation Tg . This region is described 
by the AFL in Boltzmann approximation. Under an in- 
crease of the doping, the APLD replaces Boltzmann ap- 
proximation as the valid approximation at intermediate 
temperature. In this way the transition from the AFL in 
Boltzmann approximation to the one given by the APLD 
is visible. This can be seen in Fig.dJ too, where an expan- 
sion around T = 800 K was considered and the difference 
of the chemical potential at a temperature T = 400 K was 
studied. The transition occurs therein when moving up- 
wards, respectively increasing the renormalised doping y, 
along the line for r\ = 3.4(eV) _1 . Since the transition 
temperatures To for the different doping values are al- 
ways lower than the expansion temperature T the APLD 
as described by Eq. (|2"Tj) remains valid. 



V. CONCLUSION 



This investigation of the origin of an AFL behaviour 
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Figure 8: Specific heat of CuRhi„ x Mg x C>2 and its 
derivative with respect to the temperature (inset) given 
by the calculation with the parameters of the density of 
states from Eq. (j56")l for values of the doping x = 2% 
(solid) , 3% (dashed) and 10% (dotted) . The arrows indi- 
cate the position of the characteristic (breakdown) tem- 
peratures for the different doping values. The curves 
show a transition between a Fermi liquid and an apparent 
Fermi liquid indicated by an inflexion point. 

started from a density of states of non-interacting fer- 
mions written as a Taylor series around a single discon- 
tinuous band edge. With this simplification Sommerfeld 
expansion and Boltzmann approximation could be formu- 
lated as expansion of the polylogarithm. While the ex- 
planation of Fermi liquid like behaviours at large temper- 
ature within the Fermi liquid theory would require large 
temperature scales, it can be expected from the latter ap- 
proximation for very low temperature scales, too, e.g. in 
the chemical potentials shown in the insets of Figs. 5 and 
7. However, the behaviour of this approximation presents 
distinctive differences like an offset in the thermopower 
within the TICR framework. Furthermore, its doping 
dependence is similar to the one of Heikes behaviour but 
does not describe such a temperature independent be- 
haviour. In addition, it states that the degeneracy factor 
in modified Hcikcs formulas has to be modified for broad 
bands by the fraction which is thermally accessible. 

An analytic anomaly, like the discontinuous band edge, 
further influences physical properties. Indeed it causes 
the Fermi liquid to break down before the temperature 
reaches the Fermi temperature. The missing region was 
made accessible by the APLD. The lowest order of this 
approximation treats the system at a given temperature 
with an averaged density of states and a renormaliscd 
doping. While it incorporates Boltzmann approximation, 
it becomes exact if the density of states can be described 
only by a discontinuity. In this case the large temperature 
limit results in the modified Heikes formula and therefore 
describes the temperature independent thermopower as 
given by Heikes formula correctly. Furthermore, an addi- 
tional AFL region could be identified by the APLD at a 
particular temperature. This phase strongly depends on 



the presence of a noticeable discontinuity. Similar to the 
lower tail of the Fermi function in Boltzmann approxi- 
mation, its upper tail combined with the missing states 
beyond the discontinuous band edge will then lead to an 
AFL behaviour, but at lower temperatures than the for- 
mer approximation. 

By the application of this framework on the doped dc- 
lafossites CuCr02 and CuRh02, one material could be 
identified which falls in the Boltzmann region and one 
which exhibits an AFL behaviour caused by a large dis- 
continuity. The data of the thermopower of these materi- 
als show very good agreement to the presented theory. In 
addition, the recently addressed formation of a spinel or 
CuO phase in the former material was indicated by the 
change of the parameters of the effective density of states. 
Furthermore, the determination of the effective charge 
density could be performed within the phcnomenological 
treatment the results of which were in good agreement 
with previous EDS measurements. From the latter ma- 
terial the effective doping values were determined, too, 
and a much larger discontinuity was extracted which co- 
incide with ASW2 calculation. Thereby it was seen that 
although larger doping values were considered than for 
the former material the thermopower was additionally en- 
hanced by the discontinuity. Additionally, this anomaly 
would result in a peak in the derivative of the electronic 
part of the specific heat of this material which might be 
experimentally tested. 
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